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Analytic and numerical approximations to a hydrodynamical model describing 
longitudinally expanding, cylindrically symmetric, finite systems are fitted to pre- 
liminary NA44 data measured in 200 AGeV central S + Pb reactions. The model 
describes the measured spectra and HBT radii of pions, kaons and protons, si- 
multaneously. The source is characterized by a central freeze-out temperature of 
T = 154 ± 8 ± 11 MeV, a "surface" temperature of T r = 107 ± 28 ± 18 MeV and 
by a well-developed transverse flow, (ut) = 0.53 ± 0.17 ± 0.11. The transverse geo- 
metrical radius and the mean freeze-out time are found to be Rq = 5.4 ±0.9 ±0.7 
fm and to = 5.1 ± 0.3 ± 0.3 fm/c, respectively. Fits to preliminary WA93 200 
AGeV S + Au and WA98 158 AGeV Pb ± Pb data dominated by pions indicate 
similar model parameters. The absolute normalization of the measured particle 
spectra together with the experimental determination of both the statistical and 
the systematic errors were needed to obtain successful fits. 



1 Introduction 

The reconstruction of the space-time picture of particle emission in high en- 
ergy heavy ion and particle physics became a focal point of current research 
interest. In high energy heavy ion physics, the space-time information on par- 
ticle production is needed to determine the volume of hot hadronic matter at 
freeze-out El while in particle physics Bose-Einstein correlations of pions from 
different decaying W mesons may result in large contribution to the systematic 
errors of W-mass measurements at LEP II El, a major problem in forthcoming 
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precision determination of W mass. 

In this paper, we attempt to reconstruct the space-time picture of high 
energy n heavy ion reactions by checking whether the hydrodynamical model 
of ref. U (referred to as the Buda-Lund hydro parameterization) is able to fit 
NA44 data on two-particle correlations and single-particle spectra at S + Pb 
200 AGeV central reactions at CERN SPS. Data from WA93 200 AGeV S + 
Au and WA98 158 AGeV Pb + Pb experiments are also used to confirm the 
reliability of the fit results. The Buda-Lund hydro parameterization u corre- 
sponds to a class of longitudinally expanding, cylindrically symmetric, finite 
systems, where the local rest density distribution, the local inverse tempera- 
ture distribution and the freeze-out proper-time distribution is characterized 
by their means and variances, respectively. A scaling longitudinal flow field 
is assumed together with a linear transverse flow profile. As a consequence of 
the cylindrical symmetry of the model, a l/^/roJ scaling of the HBT radius 
parameters was predicted in certain limited regions of the parameter spaceSa. 
Such a scaling is difficult to obtain in other type of models, however the scal- 
ing law is satisfied by not only the NA44 data on S+Pb reactions El but also 
by the NA22 data on hadron-proton interactions at CERN SPSH and, to some 
extent—by. preliminary e + + e~ annihilation data in two-jet events at LEP-I 0. 
In ref.Btm it was observed for the first time, that the source parameters can be 
determined precisely only if a simultaneous analysis of the particle spectra and 
correlation functions is performed. Such results on data fitting to analytical 
and numerical approximations are presented in the body of this paper. 



2 The model and its re-parameterization 

2.1 The model 

The Buda-Lund hydroparameterizationi makes difference between the central 
(core) and the outskirts (halo) regions of the collided matter. The follow- 
ing emission function S c (x,p) applies to a hydrodynamically evolving three 
dimcnsionally expanding, cylindrically symmetric finite system: 



a , w4 9 d i &(x)p IJ , 

SJX.p) d X = -TT -. - r , (1) 

i^f ( u»{x)p» fi(x)\ ' 
P \ T[x) T(x)) 

where the subscript c refers to the core, the factor d 4 T, tJ, (x)p fJ , describes the 
flux of particles through a finite, narrow layer of freeze-out hypersurfaces. 
It is assumed that any of the layers can be labelled by a unique value of 
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t = \/t 2 — z 2 , and the random variable r is characterized by a probability 
distribution, such that 

d i W(x)p l _ l = m t cosh[?7 - y] H{r)dT T a drj dr x dr y , (2) 



Here m t — ym 2 + p 2 + p 2 stands for the transverse mass, the rapidity y and 

the space-time rapidity 77 are defined as y = 0.5 log [{E + p z )/(E — p z )\ and 77 = 
0.5 log [(i + z)/(t — z)\ and the duration of particle emission is characterized by 
H(t) cx exp(— (r — to) 2 /(2At 2 )). Here To is the mean emission time, At is 
the duration of the emission in (proper) time. The four-velocity and the local 
temperature and density profile of the expanding matter is given by 

u^(x) — ( cosh [rj\ cosh [rjt], sinh[?7 t ] — , sinh^t] — , sinh[?7] cosh[rjj] ) , (3) 
V n r t I 



sinh[r? t ] = b—, r t = Jr 2 x + r 2 (4) 
T V y 

The inverse temperature profile is characterized by the central value and its 
variance in transverse and temporal direction, and we assume a Gaussian shape 
of the local density distribution: 

Ji \ ( l~ _ \2\ 

(5) 




T(x) T 2R 2 G - " ' (,! 

where n[x) is the chemical potential and T(x) is the local temperature charac- 
terizing the particle emission. Note that the strength of the transverse changes 
of the temperature profile, the gradient of the transverse flow and the strength 
of the temporal changes of the temperature profile are controlled by the di- 
mensionless parameters a, b and d, respectively. 

2.2 Analytic approximations 

In Ref. i, the Boltzmann approximation to the above emission function was 
evaluated in an analytical manner, applying approximations around the sad- 
dle point of the emission function. The resulting formulas express the Invari- 
ant Momentum Distribution (IMD) and the Bose-Einstein correlation function 
(BECElin an analytic way. The resulting analytical formulas are given in 
refs. ButrtHl and shall not be recapitulated herewith. This Boltzmann approxi- 
mation was also applied to the numerical approximate evaluation of the model, 
as given in the next section. 
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2.3 Core/halo correction 

The effective intercept parameter A*(y, mj) of the Bose-Einstein correlation 
function controls the cop.xatio in the particle production in the core/halo 
picture developed in refs.OJO'EJ. With this factor the total invariant spectrum 
in y rapidity and transverse mass mj follows as 



d n 1 d n c tt 



dy dm 2 J~\l dy dm 2 */~\ 



S c (x,p)d 4 x, (7) 



The momentum dependence of A» parameter was mesured by NA44 in ref. 
although with a very limited momentum resolution. 



2.4 Re-parameterization 

We introduce the surface temperature T r — T(r x = r y = Rq, t = tq) and the 
temperature after most of the particles were emitted as T t — T(r x — r y — 
0; r = To + y/2Ar). Here Rq stands for the transverse geometrical radius of 
the source, tq denotes the mean freeze-out time, At is the duration of the 
particle emission and we denote the temperature field by T(x). The central 
temperature at mean freeze-out time is denoted by To = T(r x = r y = 0; r = 
to). 

Then the relative transverse and temporal temperature decrease can be 
introduced as 

AT _ T - T r AT _ Tq - T t 
[-jT/r — Tj; , \-jTlt - — (o) 

and it is worthwhile to introduce the mean transverse flow as the transverse 
flow at the geometrical radius as 

K>=&— . (9) 

Hence the 3 dimensionless parameters can be re-expressed with the physical 
parameters introduced above as 

a 2 = 4- < — ) r , b = (U t ), d 2 = ^lr ( — ) t . (10) 

R 2 G x T ' ' R G x h At 2 T K ' 
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Table 1: Parameters from simultaneous fitting of preliminary particle spectra and HBT 
radius parameters with analytic and numeric approximations to the hydrodynamical core 
model. The table shows the combined results of the two sorts of fits for the NA44 S + Pb 
and the analytic fit results for the WA93 S + Au and the WA98 Pb + Pb reactions. The 
errors on WA93 and WA98 parameters are preliminary, see Section 3. 





NA44 S+Pb 


WA93 S+Au 


WA98 Pb+Pb 


Parameter 


Value 


Errors (Stat & Syst.) 


Value 


Errors 


Value 


Errors 


T [MeV] 


154 


± 8 ± 11 


154 


± 8 


146 


± 3 


to [fm/c] 


5.1 


± 0.3 ± 0.3 


4.7 


± 0.5 


5.0 


± 0.2 


R G [fm] 


5.4 


± 0.9 ± 0.7 


4.2 


± 0.5 


7.1 


± 0.3 


Ar) 


1.6 


± 0.3 ± 0.3 


1.6 


± 0.7 


1.6 


± 0.3 


At [fm/c] 


0.3 


± 0.3 ± 0.3 


0.9 


± 0.8 


1.7 


± 0.1 


a 


0.63 


± 0.09 ± 0.01 


0.21 


± 0.11 


0.06 


± 0.06 


b 


0.50 


± 0.06 ± 0.06 


0.73 


± 0.18 


0.36 


± 0.04 


d 


4.9 


± 1.8 ± 1.1 


4.7 


± 4.3 


8.1 


± 0.6 



Table 2: Parameters calculated from the results in Table 1. Preliminary fits for NA44 using 
directly these parameters indicate similar values with smaller errors. 



Calculated parameters 


Value ± Error 


Value ± Error 


Value ± Error 


m* 

{—}* 
T r [MeV] 
T t [MeV] 
<«t> 


0.44 ± 0.45 ± 0.22 
0.08 ± 0.62 ± 0.56 
107 ± 28 ± 18 
143 ± 61 ± 56 
0.53 ± 0.17 ± 0.11 


0.04 ± 0.09 
0.81 ± 13.30 
148 ± 19 
85 ± 77 
0.66 ± 0.34 


0.01 ± 0.04 
7.50 ± 3.14 
145 ± 8 
17 ± 10 
0.51 ± 0.10 



3 Fitting the model to data 

The model has been tested in two different ways using the analytical approx- 
imation referenced above and the numerical approximation that is based on 
a numerical integration of the emission function. The kinematic parameters 
are fitted simultaneously to preliminary IMD and HBT radii measured by the 
CERN NA44 experiment in central S + Pb collisions at 200 AGeV. 

As emphasized previously, absolutely normalized data were utilized in 
these fits, which were performed with the help of the CERN function mini- 
mization package MINUIT. Moreover, core/halo correction oc l/-\/A7 is applied 
and the corresponding errors are propagated properly. Due to the these con- 
ditions a unique minimum is found, and the strongly coupled, normalization 
sensitive d and At parameters are determined. In contrast, if the data are 
fitted without absolute normalization, we reproduce the results in ref. Ej and 
we obtain big errors on these parameters. 

In the numeric approximation scheme we evaluate the means and the vari- 
ances of the core as suggested in ref. EH. Since this numeric approximation 
scheme is not an exact calculation, but an approximation in a different way 
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than the analytic approach, we use it to estimate the systematic errors of the 
model parameters. 

The analytic and numeric results are combined in Table 1 to estimate the 
model parameters and their errors properly. On Figure 1, the analytic and 
numeric fits to measured data are shown simultaneously. The parameters a, 
b and d are transformed to the corresponding relative transverse temperature 
decrease, the mean transverse flow and the relative temporal temperature de- 
crease on Table 2. 

A comparision between the numerical and the analytical approximation 
schemes indicates that the minima found by the two rather different fitting 
methods coincide within 2 standard deviations for each parameter of the model. 
To estimate the systematic errors half of the difference between the minima of 
the two fits is evaluated for each parameter. To estimate the best values of the 
fit parameters, the mean of the analytic and numerical minima is taken. To 
estimate the errors on these values, the statistical error equals to the bigger of 
statistical error of the numerical and the analytical fits, the systematical error 
is defined as above. 

Fits to preliminary data of the WA93 and WA98 experiments seem to 
provide similar source parameter values like those obtained by NA44. However, 
some of the parameters have big errors because the particles are unidentified 
in these experiments and statistics allowed for determination of only one HBT 
radius in each (side, out, long) directions in the Bertsch-Pratt frame. The 
normalization of the WA98 spectrum was fixed manually, which resulted in 
artificially small errors in Table 1 . As an indication, the following characteristic 
values are obtained from fits using the analytic model approximation. In case 
of WA93 we get a central freeze-out temperature of To = 154 ± 8 MeV and a 
mean transverse flow of (u t ) — 0.66 ± 0.34. Analysis of WA98 data shows that 
the corresponding parameters are 146 ± 3 MeV and 0.51 ± 0.10, respectively. 

4 Discussion 

4-1 The Source of particles in space-time 

On Figures 2 and 3 we indicate the reconstructed space-time distribution of the 
source of particles in 200 AGeV central S + Pb reactions, as a function of the 
time variable t and the coordinate along the beam direction, z, both measured 
in the mid-rapidity frame that moves in the laboratory with yo = 3.0. The 
momentum-integrated emission function at r x = r y = is given by 

S c (t, z) oc exp {-^-r- ) oxp [-^^ ) (11) 
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Figure 1: Simultaneous fitting of particle spectra and HBT radius parameters. Full line 
shows the fit using the analytical approximation, dashed line shows the result of the numerical 
approximation. Data points on the spectra represent the jy-averaged NA44 2 dimensional 

,2 - 

2-2 — — distribution. HBT data are found in ref. . 

mt drrit dy 



where the parameters are taken from Table 1, corresponding to the best fit to 
NA44 S+Pb data in this picture. Note the relation (t, z) — (r cosh(?7), r sinh(?7)) 
The contour-plot of S c (t, z) on Figure 2 shows that the with of particle emis- 
sion is found to be narrow in proper-time. The parametric plot of S c (t, z) on 
Figure 3 indicates the long tails of particle production in the \z\ rs t regions, 
the height of the curve is proportional to the production probability. 

Note that in our case the hypothesis that pions, kaons and protons are 
emitted from the same hydrodynamical source is in a good agreement with the 
fitted data. The hypothesis that pion and kaons are emitted from different 
sources was investigated in ref. L3, and this hypothesis resulted in a worse de- 
scription of the data than the hypothesis that the pions and kaons are produced 
from the same source. 
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S(t,z) for S + Pb 
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Figure 2: Contour plot of the reconstructed source function for S+Pb reactions as measured 
by the NA44 collaboration. Note the large value of the mean freeze-out time and the narrow 
width of the emission zone in proper-time. 
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Figure 3: Parametric plot of the reconstructed source function for S + Pb reactions as 
measured by the NA44 collaboration. The source is the same as that of Figure 2, but a 
different view-point is chosen to show the shape of emission probabilities for large values of 
\z\ in the mid-rapidity frame. 



4-2 Comparision to hadron - proton reactions at CERN SPS 

The NA22 collaboration fitted the same model to 250 GeV meson + p datai. 
Their best fit parameter values were Arj = 1.36 ± 0.02, T = 140 ± 3 MeV, 
(ut) = 0.20 ± 0.07 and (££) r = 0.71 ± 0.14. From a combined HBT and 
spectrum analysis the NA22 experiment finds a mean freeze-out time of To = 
1.4 ± 0.1 fm/c, and a comparable duration parameter At > At, = 1.3 ± 0.3 
fm/c. The transverse geometrical radius was found to be Rq — 0.88±0.13/m, 
which is slightly larger than the corresponding Gaussian radius parameter of 
the proton. 

When comparing the source parameters of S + Pb reactions to h + p we 
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find that the freeze-out temperature at the mean emission time at the center 
is similar in both cases. The surface temperatures are also equal within errors 
(82 ± 7 MeV vs. 107 ± 28 ± 18 MeV). The width of the particle emission in 
spacetime-rapidity r\ is also similar, 1.36 ± 0.02 vs. 1.6 ± 0.3 ± 0.3. However, 
we find a significantly larger mean freeze-out time in S + Pb reactions as 
compared to h + p reactions (5.1 ± 0.3 ± 0.3 fm/c vs. 1.4 ± 0.1 fm/c). It 
seems that we observe a sharper freeze-out hypersurface in heavy ion collisions 
than in hadron-proton reactions, cf. At = 0.3 ± 0.3 ± 0.3 vs. At > 1.3 ± 0.3 
fm/c. The transversal radius of the matter is also larger in S + Pb reactions, 
5.4 ± 0.9 ± 0.7 vs. 0.88 ± 0.13 fm. We also find a much stronger transverse flow 
in S + Pb than h + p reactions. 



4-3 Comparision to other NA44 data 

A comparisicp. of this analysis of NA44 data for S + Pb reactions to recent 
NA44 resultsllll on p + p, S + S and Pb + Pb reactions shows that the freeze- 
out temperature in the center at the mean freeze-out time, T , in our case of 
S+Pb reactions is within errors similar to the values obtainedJjrom an analysis 
of the particle spectra of p + p, S + S and Pb + Pb reactions O. In the S + Pb 
case, we have a complete description, including the low pt part of the pion 
spectra with an acceptible x 1 /ndf due to our inclusion of the temperature 
inhomogeneities. Note also that we do not simply reproduce the slopes of 
particle spectra but we reproduce the absolute normalized particle spectra for 
pions, kaons, as well as the full shape of the proton spectra and the HBT radii. 



5 Conclusions 

A combined fitting of HBT radius parameters and spectrum fitting is presented 
in the Buda-Lund hydrodynamical parameterization scheme for CERN SPS 
heavy ion reactions. 

It was necessary to use different particles and absolutely normalized par- 
ticle spectra to obtain a well-defined minimum of each parameters for data 
from the small NA44 acceptance. The model is shown to describe the data 
in a statistically acceptable manner. The parameters can be determined with 
an accuracy of 10 - 20 % relative errors including systematic uncertainties. 
The transverse flow and the radial temperature inhomogeneity is necessary to 
achieve this result. 

The same model describes meson - p spectra and correlations at CERN SPS 
with a similar central temperature, similar surface temperature and similar 
space-time rapidity width. The main difference between S + Pb and (ir/K) +p 
reactions is that the geometrical radii, the freeze-out time and the transverse 
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flow are much larger in the first case. The particle emission seems to happen 
more suddenly in S + Pb as compared to hadron induced reactions at the same 
energy range. 
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